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A Laplacian Expansion for Hermitian-Laplace Functions of 

High Order* 

By E. C. MOLINA 

Among the wide variety of practical and theoretical problems con- 
fronting the telephone engineer, there is a surprisingly large number to 
whose solution mathematics has made notable contribution. In his kit of 
mathematical tools the theory of probability is a frequently used and most 
effective instrument. This theory of probability contains a large number 
of theorems, a large number of functions, which permit of application to 
telephony. Among these is a particular tool, a particular group of mathe- 
matical functions known as the "Hermitian Functions," each of which is 
identified by a number called its "order." These mathematical functions 
or relations have no practical utility until the variables in the equation can 
be assigned numerical values and the resultant numerical value of the 
function calculated. Tables of the numerical values of Hermitian functions 
of low order exist; for example, Glover's Tables of Applied Mathematics 
cover the ground for those of the first eight orders. But tables for the 
functions of higher order are still a desideratum. This paper presents an 
expansion by means of which the evaluation of a high order function can 
be readily accomplished with a considerable degree of accuracy. 

The development of the expansion is prefaced by some remarks on the 
early history of the Hermitian functions and the relation of this history 
to modern theoretical physics. 

I 

AMONG contributions made by Laplace to the domain of pure and 
applied mathematics, two of great practical value are; 

(a) His method of evaluating definite integrals f whose integrands 

involve factors raised to high powers ; 

(b) The pair of orthogonal polynomial functions 2 which he defined by 

the following Equations (1) and (2) 

(1) [(2»)!V^/2 2 »»!]Z7„(«) = r e-' 2 (x - iu) 2 "dx 

•7 — 00 

f*QO 

= 2e" 2 I g-*V B cos (2ux)dx; 
Jo 

e~*(x - iuy^dx 

00 

Jt*a> 
e -*V n+1 sin (2ux)dx. 
o 

* Presented at International Congress of Mathematicians, Oslo, Norway, July 
13-18, 1936. 
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These polynomial functions formed the coefficients in a series sat- 
isfying the partial differential equation 

-, = 2U+2u- + -^ 

to which Laplace reduced the solution of the following ball problem : 3 

Consider two urns A and B each containing n balls and suppose that of the total 
number of balls, 2n, as many are white as black. Conceive that we draw simultane- 
ously a ball from each of the urns, and that then we place in each urn the ball drawn 
from the other. Suppose that we repeat this operation any number, r, of times, 
each time shaking well the urns in order that the balls be thoroughly mixed: and 
let us find the probability that after the r operations the number of white balls in 
urn A be x. 

Under the caption "The Statistical Meaning of Irreversibility" 
Lotka 4 has pointed out the significance of Laplace's ball problem in 
the modern kinetic theory of matter. Moreover, Hostinsky B has 
shown the bearing of the same problem on the theory of Brownian 
movements and said "In effect, the partial differential equation ob- 
tained by Laplace has been refound by Smoluchowski." 

To avoid confusion with the Laplace functions which one encounters 
in spherical harmonic analysis, the functions defined by Equations (1) 
and (2) are herein designated as Hermitian-Laplace functions. Such 
a designation is justified by the Equations (3) and (4) derived in the 
next paragraph. 

II 

We also find in Laplace 6 

/»(«)* - f V*V" cos (2ux)dx = ( "ff T (^2^) . 

r°° (— l) n+1 ->fTr (d 2n+1 e~ u2 ) 

/«'(«) = I *-*V n+1 Sin (.2liX)dx = 2 2n + 2 du 2n + i • 

Jo 

Comparing these Laplacian expressions for the definite integrals /„(«) 
and /»'(») with the Equations (1) and (2) we see immediately that 

(3) £/„(«) = (- l)*[»!/(2»)!]fltn(«), 

(4) U n '(u) = (- l)«+ 1 [«!/2(2» + l)!]tf 2n+ iO), 

where H 2n and H 2n +i are the original Hermite polynomials 7 of order 
In and 2« + 1, respectively. These equations connecting the Her- 

* The symbols I n (u) and /„'(«) are introduced here as convenient abbreviations 
for the integrals to which they are equated; these symbols do not appear in Laplace. 
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mite with the Laplace polynomials have been presented in an earlier 
paper. 8 

Appell and Feriet, Arne Fisher, T. C. Fry, H. L. Rietz and others 
base their definitions of the Hermite polynomials on e - * 2 ' 2 instead of 
e~ xi . We shall write A n (u) for the wth polynomial as defined by these 
authors, reserving H n (u) to symbolize the Hermitian polynomial as 
defined in his paper of 1864. Thus, in what follows, 

A n (x) = (- \) n e x<i i 2 (d n e- x2 i 2 Jdx n ), H n (u) = (- lV2)M n («V~ 

III 

Laplacian expansions* for the U, H, and A polynomials follow 
immediately from those obtainable by applying to the integrals /„(«) 
and I„'(u) his method of evaluating definite integrals whose integrands 
embrace factors raised to high powers. As will be shown in Part IV 
of this paper, we have 

I n (u)/[_^(Y^n = [Scos {uj2N)l + [5' sin (W2JV], N = 2», 
/»'(«)/[ V^(FViV) w ] = IS sin (mV2N)] - IS' cos (mV2W], N - In + 1, 
where Y = (xe~' 2 ) for x = X = 1/V2 and 



The explicit expressions for K , K 2 , K A and K h K a , K b are given in 
Section V of this paper. The desired expansions are then given by 
the equations: 

e» 2 I„(u)/^r = f/ n («)[(2«)!/2 2 »+ 1 »!] 

= H 2n {u)\_{- l)»/2 2 «+i] 

= A tn (uj2)l(- l)"/2 n+1 ], 
e» 2 I n '(u)/^= U n \u) [(2« + l)!^ 2 -*^!] 

- 7/ 2n+1 («)[(- 1)»+V2 2 " +2 ] 

= ^ 2n+1 («V2)[(- l)"/2" +, V2]. 

The numerical results shown below in Table I indicate the efficacy 

* It may be of interest to compare the expansions presented in this paper with 
the asymptotic forms of the Hermite functions given by N. Schwid 9 and by M. 
Plancherel and M. Rotach. 10 
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TABLE I 



N 




True Value of Aj>j(x) 


l,A N (x) - A N (,x)]/A N (x) 






s = 1 


s = 2 


s = 3 


9 


0.1 
0.5 
1 

2 
3 
4 


9.32438 X 10 1 
3.26533 X 10* 
2.80000 X 10 1 

- 1.90000 X 10 2 
1.62000 X 10 3 

- 1.74680 X 10 4 


0.0089 
0.0084 
0.0263 
0.0002 

- 0.0474 

- 0.0283 


0.0000 
0.0001 

- 0.0010 
0.0018 

- 0.0027 

- 0.0082 


0.0000 
0.0000 

-0.0002 
0.0001 
0.0002 

- 0.0027 


10 


0.1 

0.5 

1 

2 

3 

4 


- 8.98064 X 10 2 
4.90439 X 10 1 
1.21600 X 10 3 

- 2.62100 X 10 3 
9.50400 X 10 3 

- 5.18090 X 10* 


0.0084 

- 0.0027 
0.0046 

- 0.0147 

- 0.0436 
0.0445 


- 0.0001 
0.0013 
0.0003 
0.0002 

- 0.0044 
0.0013 


0.0000 
0.0000 
0.0000 
0.0001 

- 0.0004 

- 0.0018 


15 


0.1 

0.5 

1 

2 

3 

4 


- 1.98001 X 10 s 

- 5.05845 X 10 5 
4.69456 X 10 5 

- 1.41980 X 10 s 
4.38955 X 10 6 

- 1.85644 X 10 7 


0.0054 
0.0052 
0.0022 

- 0.0102 

- 0.0284 

- 0.0338 


0.0000 
0.0000 
0.0002 
0.0000 

- 0.0017 

- 0.0041 


0.0000 
0.0000 
0.0000 
0.0000 

- 0.0001 

- 0.0006 


20 


n.i 

0.5 

1 

2 
3 

4 


5.90233 X 10 8 

- 4.45178 X 10 8 

- 1.61935 X 10 8 

- 1.62882 X 10 9 
4.60718 X 10 9 
8.53219 X 10 9 


0.0042 
0.0035 
0.0046" 
- 0.0081 
-0.0212 
0.1241 


0.0000 
0.0000 
0.0000 
0.0000 
- 0.0009 
0.0068 


0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 



of the Laplacian expansion as applied to the evaluation of A n (x), 
x = ma/2, for values of x ranging from 0.1 to 4 and for N = 9, 15, 10 
and 20, respectively. 

Designating by 8 An{x) the approximate value obtained for An(x) 
when one takes into account the first 5 terms in each of the two series 
S, S', the last three columns of the table show the proportional errors 
incurred when 5 = 1,2 and 3, respectively. It will be noted that for 
N = 10 and x = 4 the second approximation is closer than the third; 
this situation will occasion no surprise if it is recalled that to obtain 
the best results the natural order of the terms may have to be altered 
when, for example, one expresses a term of the binomial expansion in 
a series of Hermite functions. 

For the convenience of one who wishes to calculate An(x) for values 
of N other than 9, 10, 15 and 20, there are given in Table II the values 
of ur< m+1) K m for m = 0, 1, 2, 3, 4, 5 and those values of u covered by 
Table I. 

I am indebted to Miss E. V. Wyckoff of Bell Telephone Labora- 
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tories, for the computations involved in the preparations of Tables I 

and II. 

TABLE II 



uV2 


m-'Ko 


u-»Ki 


u-«K« 


0.1 

0.5 

1 

2 

3 

4 


0.7053412 

0.6642654 

0.5506953 

0.2601300 

0.07452849 

0.01295111 


0.234229 
0.198970 
0.0936947 

- 0.158968 

- 0.121885 
0.0244632 


0.034490 

- 0.069716 

- 0.273541 

- 0.018172 
0.514828 

- 0.078479 




u~*Ki 


u~*K, 


u-»Ki 


0.1 
0.5 

1 
2 
3 
4 


- 0.03520828 

- 0.1591469 

- 0.2294564 

- 0.08671002 
0.05589638 
0.04317037 


0.00602953 
0.0450683 
0.150921 
0.255634 

- 0.104689 

- 0.162097 


0.091905 

0.40870 

0.48658 

- 0.78050 

- 0.35023 
0.98512 



IV 



A simple change of variable gives 

/»00 

(5) /„(«) = (^) N+1 (e-*x) y cos (x2u^N)dx, N = 2«, 

(6) !„'(«) = (4N) N+l (r*x)*mn (x2ujN)dx, N = 2» + 1. 

Jo 

Set y(x) = e~ x *x, and note that dy/dx = for x = X = I/V2T Now set 

Y = y(X), 

Lg(x)J = (log Y - \ogy)/(x - Xf 

--w-im+im-im +■■■]■ 

(7), t = (x - X)g(x). 
These transformations give 

/•CO 

/„(«) = (-y[N) N+l Y N I e~ Nl \dx/dt) cos (xlujN)dt. 

1/ — * 

By (7) and the Lagrange-Laplace expansion for a function of x in 
powers of t we obtain 

!„(«) = tfN)» +1 Y« t r e- Nl V 2m ^2j(2m)Udt 

m=0 J -co 
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or 

I n (u)l^{Y4N) N = £ (l/2V^) 2m (^2 m /w!), 

m=0 

where, writing D for the differential operator d/dx, 

A 2m = lD*»r lM) cos {xlujN)^ 
or, by the Leibnitz theorem for the product of two functions, 

A* m = E o ( 2 ^)(2«ViV) r cos (u^2N + nr/2)lD x *°~r itm+1) l*-x 
and, therefore, " 

^ = cos (um) g(% y- m «W LD.~r~»^ 

- sin (WSV) g^ ^-'^y^' CW— «-—>!- 

on separating the even and odd terms in r. Now setting m — r = s 
and summing with reference to 5 and r, instead of m and r, gives 



£ (l/2ViV) 2 ^ 2m /m! 

l_y. ! « (2r + 2s)! (- w 2 )' 



m=0 



/> OT\^/ * V* 1 ^ (2r + 2s)! (-u 2 Y 

x ^r;™]^ _ sin (W 2iv) £ (^)^_L_ 

A rt-o (2r + l)l (r + j)! L ^* g J *= x " 

But, writing w/g = v, we have 

r oo (_ 1 •\r+»„2r+2»"| 

= (_ 11««-(2«+l) D 2« 7 ,2«+in 2* V ^ ' 

= (- l)'M-»«+ 1 >[ZV«» 2H - 1 ZV«e- , ' a ] I 



since D v 2, v 2m = for m less than 5. 
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Likewise 



« (2r + 2*)l (- 1)^+' VD w (wrfD-i 
Therefore, finally, 



/.(«) 



/ l y 



(wt)% - cos (W2 "> £W) [**~»*J<- D« 



where 



+ sin (»m(^=)± •fc$™ c«- 2 ^-](- 1)- 1 . 

(25) IX*. = [Z)x 2 V' +1 e-" 2 i72.(y)]x=x, 
(25 - 1) !£*._i = [Z?, 2 -V«- a Jr,._i(w)]^x. 

Substituting sin (xlu^N) for cos (x2u4N) in the equations defining 
.4 2 m and then proceeding exactly as above we derive the corresponding 
expansion for I n '(u). 

V 
To obtain the values of K u and -K" 2s _i note that 

Xg* = 1 - (* - X)/3X + (x - X) 2 /4X 2 - (x - X) 3 /5X 3 + • • • 
gives, for x = X = 1/V2, 

v = (1/V2)«, 
dw/dx = (l/6)w, 

</y</* 2 = - (1/3V2)«, 
<ft/4* = (53/90)m, 
<*y<fcc 4 = - (21lV2/13S)«, 
dVrfx 5 = (79/7)«, 



J2(u- l K Q ) = <r»« 2 , 

36J2(u-*K 2 ) = e-*» 2 (tt 6 - 6m 4 - 9m 2 + 12), 

7776V2(«- 5 i^4) = <H U V 2 - 12m 10 - 183.6m 8 + 1432.8m 6 

+ 2889m 4 - 10368m 2 + 432), etc. 

6(m- 2 .Ki) = iur**\u % - 3), 





S = 


V2, 




dg/dx = 


- 1/3, 




d 2 gfdx 2 = 


(4V2)/9, 




d 3 g/dx* = 


-88/45, 




d*g/dx* = 


(824V2)/135, 




dfg/dx 6 = 


- 28184/567, 


etc. 






Therefore 
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648(m- 4 ^T 3 ) = we-*" 2 (M 8 - 9m 6 - 59.4a 4 + 279m 2 + 54), 

233280(m- 6 A: 5 ) =wt*\u u - 15m 12 - 414m 10 + 4494m 8 

+ 25152.4m 6 - 168723m 4 - 119340m 2 + 304560). 
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